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Abstract. A review of the known different methods and results devised to 
study the two-electron atom problem, appeared in the early years of quantum 
mechanics, is given, with particular reference to the calculations of the ground 
state energy of helium. This is supplemented by several, unpublished results 
obtained around the same years by Ettorc Majorana, which results did not con- 
vey in his published papers on the argument, and thus remained unknown until 
now. Particularly interesting, even for current research in atomic and nuclear 
physics, is a general variant of the variational method, developed by Majorana 
in order to take directly into account, already in the trial wavcfunction, the 
action of the full Hamiltonian operator of a given quantum system. Moreover, 
notable calculations specialized to the study of the two-electron problem show 
the introduction of the remarkable concept of an effective nuclear charge dif- 
ferent for the two electrons (thus generalizing previous known results), and 
an application of the perturbative method, where the atomic number Z was 
treated effectively as a continuous variable, contributions to the ground state 
energy of an atom with given Z coming also from any other Z. Instead, con- 
tributions relevant mainly for pedagogical reasons count simple broad range 
estimates of the helium ionization potential, obtained by suitable choices for 
the wavefunction, as well as a simple alternative to Hylleraas' method, which 
led Majorana to first order calculations comparable in accuracy with well- 
known order 11 results derived, in turn, by Hylleraas. 



1. Introduction 

It is commonly believed that the early wide acceptance of quantum mechanics was 
triggered by the successful description of the simplest atomic system, i.e. the hydro- 
gen atom, according to already known spectroscopic data. This is usually deduced 
from the fact that the analytic solution for the wavefunction of the (isolated) hy- 
drogen atom indeed played a pivotal role in subsequent applications of the new 
quantum theory, during the early part of the XX century. However, historians of 
physics know quite well [1] that the definitive abandonment of the old quantum 
theory by N. Bohr and A. Sommcrfeld came through the thorough consideration of 
the helium atom problem, with particular reference to the specific prediction of the 
ground state energy of neutral helium (or other measured quantities related to it). 

In 1913, following the discovery of the atomic nucleus [2], Bohr succeeded [3] in 
explaining the energy levels of the hydrogen atom in terms of quantization of the 
action for the classical Kepler orbits. Numerous attempts to explain the ground 
state of helium bu quantizing different two-electron periodic orbits in a similar 
manner, however, failed. For example, Bohr firstly discussed a simple model where 
both electrons in the helium atom move along the same circular orbit and are 
located in the opposite ends of a diameter [3]. This followed the 1904 proof by the 
Japanese H. Nagaoka [4] (obviously in the framework of classical mechanics) that 
such motion is mechanically stable (for sufficiently large attractive forces), and does 
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have the lowest possible energy. In general, the attempts to quantize the helium 
atom in the Bohr-Sommerfeld theory [5] were always based on the assumptions that 
the ground state is related to a single periodic orbit of the electron pair, and that 
the electrons move on symmetric orbits with equal radii for all times. Furthermore, 
it was as well assumed that orbits where the electrons hit the nucleus are excluded, 
and that the quantum number n in the quantization condition is an integer. Such 
approach proved unsuccessful mainly because it failed to explain the diamagnctism 
of a two-electron atom in its lowest energy state and the value of the obtained 
energy was only in a rather poor agreement with experiments. 

As reported by Briggs [BJ, in a letter to Sommerfeld of 1922, W. Heisenberg 
sketched a model where the electrons moved asymmetrically in opposite directions, 
and the quantum number associated to the quantization of the relative motion of 
the orbits under the influence of the mutual electron interaction was allowed to be 
half-integer. Interestingly enough, Heisenberg predicted a value for the measured 
ionization potential very close to that claimed by experiments of the time [7j , but 
did not publish his results due to the criticism of Bohr and W. Pauli about the 
concept of non-integer quantum numbers. However, it was quite soon recognized 
that Heisenberg' success was accidental: similar attempts to describe excited states 
of the helium [5], indeed, failed. It then became clear that "some radical modifica- 
tion in the ordinary conceptions of the quantum theory or of the electron may be 
necessary" [5J. 

The failure of the old quantum theory to describe successfully two-electron 
atoms, in fact, triggered (at least in part) the development of quantum mechan- 
ics in 1920s [T] and, once the basic formalism had been established by Heisenberg 
and E. Schrodingcr, early variational calculations (the problem couldn't be solved 
analytically, as instead for the hydrogen atom) produced remarkably good results 
for the ground state of the helium atom, thus breaking the ground for the general 
acceptance of quantum mechanics. 

In the present paper we review the known early results about this topic (Sect. 
[5]), and then discuss unknown results obtained almost simultaneously by the Italian 
physicist Ettore Majorana (Sect. [3J. As we will see, a large part of Majorana's 
results were deduced by making recourse to novel methods not yet appeared in the 
literature (both of that time and in present day studies), and while part of those 
numerical results were (and are) inaccurate when compared with the experimental 
data, the novel methods reveal quite useful in the frontier research related to atomic 
and nuclear physics. Finally, in Sect. 3] we give our conclusions and outlook. 

2. Solutions to the helium atom problem 

Within the Schrodingcr framework of quantum mechanics, atoms with two electrons 
are described by a wavefunction ip satisfying the eigenvalue equation (in electronic 
units): 

i?V = -V 2 V + 2| -—- — + —) ^ = W^, (1) 
V r i r 2 r 12 J 

where n, are the distances of the first and second electron from the nucleus, 7*12 
their mutual separation and W the energy of the system in Rydberg; V 2 = V 2 + V 2 ., 
with V 2 , V 2 , the Laplacian operators in the space of the first and second electron, 
respectively This differential equation is not separable, so that, unlike the solutions 
for the hydrogen atom, its solutions for the eigenfunctions and energy eigenvalues 
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cannot be expressed in closed analytic form, and thus approximation methods must 
be used. As it became clear soon, the key point where different strategies could 
be tested (thus proving or disproving the underlying theory) was the prediction 
about the ground state energy - mainly for the neutral helium atom - so that 
in the following we focus on this issue. Quite an exhaustive, though apparently 
outdated, review of the complete two-electron problem may be found in Ref. |10j . 
For future reference, it is useful to recall the experimental value of the ionization 
energy Wi = —Z 2 — W (here it is sufficient the accuracy as in late 1920s - early 
1930s): 

Wj xp = 24.46 cV . (2) 

2.1. Perturbative calculations. The most simple method to solve Eq. ([I} is to 

find a soluble problem to serve as an unperturbed solution to which the theory of 
perturbations may be applied. This was carried on by Sommerfeld's student A. 
Unsold in 1926 [11], by factorizing the wavefunction in terms of spherical coordi- 
nates. The main idea and the basic calculations may be more easily understood 
by following a bit later derivation by Majorana, now reported in Ref. [12]. In first 
approximation, the interaction among the two electrons may be neglected, so that 
the eigenfunction of the unperturbed hamiltonian, 

H = -V^-2(- + -) , (3) 

may be expressed as a product of two hydrogenic Is states, 

i> = c- Zri c- Zr2 , (4) 

with energy W = —2Z 2 . The inter-electron potential is, then, treated as a pertur- 
bation: by applying standard computational techniques, the first order result for 
the ground state energy (or ionization energy) is the following: 



W = -2Z 2 + -Z 
4 



9 5 

Wi = Z 2 - -Z 



(5) 



For the neutral helium and ionized lithium, this translate in the predictions W^ e = 
20.3 eV and W\ K = 71.0 cV, respectively [Tl]Q As noted by J.C. Slater few months 
after the appearance of the Unsold results [13], "the solutions so far obtained, 
though qualitatively interesting, are quantitatively far from the truth" when com- 
pared (for example, for the helium) with the experimental result This urged 
the same Slater to search for an alternative strategy to the problem, by changing 
the starting point (that is, the unperturbed system) of the perturbation theory. 
The motivation is, apparently, mathematical in nature: 

The method so far used in such cases has been to leave out, for 
the unperturbed problem, all those terms in the energy which pre- 
vent the separation of variables, regarding all these as perturbative 
terms, and solve the remaining problem by separation. Since these 
are, in general, the terms representing the interaction between elec- 
trons, it is obvious that the unperturbed problem will be a very poor 
approximation to the real case. We proceed in quite a different way, 
making no attempt to separate variables at all j!3j . 



^Or, more precisely, according to Majorana, W^ e = 20.31 eV and W\ l = 71.08 eV. 
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The idea was to treat the variables of one of the two electrons as parameters in 
the Schrodinger equation of the other electron (that is, the equation for an electron 
in the field of the rest of the atom). From a more physical perspective, Slater 
imagined one electron in a fixed position and then obtained for the second electron 
a two-center problem, the eigenfunctions of which could be evaluated numerically as 
function of the distance between the two centers. For the first electron he deduced, 
thus, a numerical equation in which the eigenvalue of the remaining electron was 
introduced and which represented the motion of this electron in a central field. 
Then, both solutions were integrated into an approximate solution for the two- 
electron problem, and finally a perturbative calculation was performed. By using, 
initially, poorly accurate solutions to the unperturbed problem, Slater found a value 
Wj 10 = 24.35 eV, apparently very close to the experimental determination, but in a 
subsequent note added in proof [13] he corrected such a prediction, by obtaining a 
numerical value of 5-^6% greater than the experimental value (which was, however, 
a better result when compared with the Unsold prediction). Almost a year later, 
Slater published a paper |14j where he analyzed the reasons of why his perturbative 
method got worse with increasing perturbative order, despite the outer electron in 
the helium atom was assumed to move not in a hydrogen-like field, as for Eq. ([4]). 
but in a certainly more accurate central field. The wavefunction replacing that in 
Eq. (|3]) (for the helium problem) for this purpose was the following: 

ip = c - 2(ri+r * )+ri * . (6) 

Slater found that his method indeed gave a good approximation to the energy levels 
and to the perturbed wavefunction, except at very small distances. More in general: 

Almost any method of setting up approximate wavefunctions will 
give fairly good results in the non-penetrating part of the orbit; 
and it is this part that is of almost complete importance in many 
practical applications that demand knowledge of the wavefunction 
[...] The operator H, almost alone among the operators of physical 
interest, puts a sever strain on the most inaccurate part of the wave- 
function, that connected with the penetrating part of the orbit, on 
account of the differentiation and the factors 1/r that it contains. 
For this reason, a wavefunction which is decidedly satisfactory for 
most physical applications may give a very inaccurate matrix of H. 
Yet it is just this matrix that is used in the perturbation method for 
solving the atomic structure problem. It seems necessary to con- 
clude from this that the perturbation method is not a satisfactory 
one for actual calculations [14]. 

A different method, then, appeared to be required for improving the predictions of 
quantum mechanics about two-electron atomsQ 

2.2. Ritz variational method. Such a novel method came into the play a bit 
earlier by means of G.W. Kellner [15], a student of M. von Laue in Berlin, who 
introduced for the first time the variational Ritz method |16] in order to estimate 
the ground state energy of the helium atom. This is given by the minimum value 



Note, however, that Slater applied his conclusions to any atom. 
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with respect to any wavefunction ip, but, when considering only a given set of func- 
tions i/>, the minimum will correspond to an approximate value, which is better 
and better as the set is enlarged. When this set reduces to the unperturbed wave- 
function in Eq. Q, the variational method gives the same result for the ground 
state energy seen above, so that a larger set of functions is required for a better 
approximation. This is achieved by choosing a trial wavefunction depending on 
some arbitrary parameters, the values of which are determined by requiring the 
energy functional in Eq. (|7|) to be the lowest possible. This minimum value then 
gives an upper limit for the ground state energy, whose approximation is much bet- 
ter as the trial wavefunction is close to the "true" eigenfunction of the problem, a 
suitable form for that being set up from given physical considerations. In principle, 
by including a sufficient number of arbitrary parameters, one can approximate the 
eigenfunction and eigenvalue of the problem as closely as desired, but it is clear 
that an appropriate choice of the form of the wavefunction will make the procedure 
converge rapidly. 

Kellner introduced the novel idea of an effective nuclear charge experienced by 
one of the two electrons in the field of the nucleus and of the other electron: this 
electron somewhat shields the nucleus, thus reducing the effective charge. This 
quantity was introduced as an arbitrary parameter in the trial wavefunction re- 
placing that in Eq. By following, again, the simpler (and somewhat more 
general) later reworking by Majorana [12] , the wavefunction thus takes the form 

.0 = c -fc(n+r 2 ) (g) 

and, minimizing the energy functional in Eq. ([7]), the following value of the k 
parameter is obtained: 

*-*-£■ (9) 

For the helium atom, then, the effective nuclear charge is less than two, kn e = 
27/16, corresponding to a ionization energy W^ e = 217/128 (in Rydbcrg units) or 
Wj 10 = 22.95 eV which is in a better agreement with experiments when compared 
with the values obtained with perturbation methods. Majorana |12j also found a 
general formula for the ground state energy (or ionization energy) by using the 
present method: 
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Wl = z 2 --z + — 

4 128 



(10) 



Instead, turning back to the helium atom, Kellner [TS] also re-iterated the method, 
by using more involved wave functions, where the exponential factor in ((5J) was 
multiplied by different combinations of Laguerre polynomials (in r\ and t<l) and 
spherical harmonics (in the cosine of the angle between n and r<z). After four 
approximation steps, the better result he obtained was Wf^ c = 23.75 eV, thus 
showing the Ritz variational method to be suitable for accurate quantum mechanical 
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treatment of the ground state of helium, depending on appropriate choices for the 
wavefunctions. 

2.3. Hartree's self-consistent field method. A different idea in managing, in 
general, many-electron atoms was introduced by D.R. Hartrce at the end of 1927 
[17j . He started by assuming the wavefunction of the given atom to be written in 
the form of a product of single particle wavefunctions y%{ri), one for each of the 
electrons, such as, for example, the one in (|4]) for the helium atom. The potential 
energy acting on the first electron is taken to be the Coulomb field of the nucleus 
corrected by the potential of the charge distribution [^(t^)) 2 of the second electron, 
averaged over all positions of the second electron: 



This corrected central potential was, then, substituted into the Schrodinger equation 
for the first electron and solved for the wavefunction yi(ri). A similar procedure is 
carried out for the second electron and, from the solutions y'^ri) for all electrons, 
Hartrce calculated a distribution of charge to be inserted again in Eq. (fTTj) to 
obtain the field of the nucleus and this charge distribution. In general, this "final" 
field is not the same as the "initial" field, since the initial, trial wavefunctions yiiu) 
do not agree with the final ones y'^Ti). The whole procedure is, then, repeated by 
using the final field of the first approximation as the initial field of the second one, 
and this is repeated over and over again until the initial and final wavefunctions 
agree to the desired accuracy. The final field is thus the same as the field produced 
by the charge distribution of the electrons, so that the field in which the electrons 
are assumed to move when calculating atomic quantities of interest (such as the 
ground state energy) is "self-consistent" . 

Here, Hartree's basic assumption was that many-electron atoms could be de- 
scribed by wavefunctions given as the simple product of one-electron wavefunctions 
of individual electrons, rather than as one complete wavefunction (as, for example, 
that in Eq. ©). This was criticized by J. A. Gaunt Q2] and others but, although 
Hartree's method "seemed to many persons to stand rather apart from the main 
current of quantum theory" since it "seemed to contain arbitrary and empirical 
elements" [19], Gaunt admitted that "Hartree's wavefunctions have been shown to 
be good approximations" [18]. This was clarified a bit later by Slater [20], who 
evaluated the error bars of Hartree's results, by making a (general, not limited 
to helium) perturbativc calculation of the energy levels by using the wavefunction 
obtained with the Hartree method. 

Hartree applied his method to Rb, Rb + , Na + , Cl~ and, first of all, to helium, 
obtaining a value for the ionization potential, Wj 10 = 24.85 eV, very close to the 
experimental determination 

In passing, we also mention a different, independent attempt made by Slater 
|21j to find an accurate expression for the wavefunction of the ground state of 
(ortho-) helium, which was a further refinement of previously considered Eq. (0). 
The general idea was to use theoretical determined functions for the limiting cases 
of large and small r-distances and then interpolate between them. In particular, 
when one electron is a considerable distance from the nucleus (large r2), while the 
other close up (small n ) , the wavefunction is assumed to be the product of an ionic 
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function of the inner electron and a hydrogen-like function of the outer one: 



(the opposite case, with small r 2 but large r±, is obtained by interchanging r\ and 
T2 in this expression). Instead, when both electrons are close up, the wavefunction 
is assumed to be of the same form as in ([6]) , but with a further correction term: 



The novel term weigthed by d was added by Slater in order to have a smooth 
interpolation between the two expressions in Eq.s (fP2"]) and (|T3|) . for intermediate 
values of r\,r%- Note that such wavefunctions were not intended for evaluating 
the energy eigenvalues, since the coefficients a, 6, c, d were just deduced from the 
experimental spectroscopic data. Instead, they were introduced in order to com- 
pute several other atomic properties of helium, such as diamagnetic susceptibility, 
repulsive forces between two helium atoms, and so on. According to Slater [2Tj : 
The calculations were made before the writer saw Hartree's pa- 
per, in which he obtains a charge density distribution for helium 
in quite a different way. The wavefunction found in the present 
paper is more complicated than Hartree's in the matter of the way 
in which it takes the interaction energy between electrons into ac- 
count. But the charge density can be computed equally well from 
either method, and this permits a comparison of the present results 
with Hartree's. The discrepancies between the two are nowhere 
greater than one or two percent. This is highly satisfactory, both 
in that it verifies the present method and Hartree's, and also that 
it justifies us in believing this density distribution to be correct 
within a narrow limit of error. The other numerical results are also 
gratifying. 

In a sense, the trial stage of the novel theory of quantum mechanics was completed, 
and the subsequent need was just to have a "standard" (and mathematically reliable 
for numerical purposes) procedure for evaluating atomic quantities to the desired 
accuracy. For many-electron atoms, this was just found in the Hartree method, as 
modified by V. Fock [22] to take exchange effects into account (fully symmetrized 
wavefunctions in the form of products). Instead, for helium, the story was somewhat 
different. 

2.4. Hylleraas variables. Early in 1928, important novelties came into play in 
the variational approach to the ground state of the helium atom from E.A. Hyller- 
aas, a Norwegian student of M. Born in Gottingen [23]. The variational method 
already exploited by Kellner was adopted because, according to Hylleraas, a "use- 
ful feature" of it was that "the eigenvalues obtained in successive orders decrease 
monotonically, and that therefore an exact eigenvalue can be given on the approach 
from one side" [23]. In the first of a series of papers devoted to the subject, he soon 
noticed [21] that the wavefunction of the two-electron problem, which, in general, 
would depend on the six coordinates ri, r 2 , actually depended only on a subset of 
three, namely ri, 7-2,0 (8 being the angle between the two vectors ri,r2), defining 
the shape of the electron-nucleus triangle. Its absolute orientation in space was, 
indeed, out of interest for the problem at hand. This introduced an enormous 
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simplification into the variational problem: indeed, Hylleraas himself was able to 
perform computations till to order 11 (instead of order 4 calculations performed 
by Kellner) , thus obtaining an excellent approximation for the ionization energy of 
helium, Wf c = 24.35 eV. 

The predictions by Kellner [15], Slater [13] and Hylleraas [24] agree within the 
error bars for direct measurements (by electron collision), but they resulted to be 
not precise enough in comparison to much more precise spectroscopic measurements 
by T. Lyman [35] ■ Hylleraas was well aware of this fact: 

The end result of my calculations [...] was greatly admired and 
thought of as almost a proof of the validity of wave mechanics, 
also, in the strict numerical sense. The truth about it, however, 
was, in fact, that its deviation from the experimental values by an 
amount of one-tenth of an electron volt was on the spectroscopic 
scale quite a substantial quantity and might as well have been taken 
to be a disproof [23] . 

For this reason, he continued to search for possible improvements in the subsequent 

years. 

In a second paper |26j on the topic, instead of adopting the r\, r 2 , 6 coordinates, 
he chose only "metric quantities appearing immediately in the expression for the 
potential energy", that is r\,r2,ri2- In particular, he introduced the notations 
(later known as "Hylleraas variables" ) 

s = n + r 2 , t = r 2 -n, u = r 12 , (14) 

so that the wavefunction had to be written just in terms of these coordinates. The 
variational method was then applied to the function -0 = (p(ks, kt, ku) - that is, 
with fc as a variational parameter - with 

<p(s,t,u) = c- s/2 P(s,t,u), (15) 

P(s, t, u) being a polynomial in s, t 1 u variables. Note that already Kellner had tried 
to force a faster convergence by multiplying the metric arguments in his wavcfunc- 
tions by a free constant, but here the novelty was in the choice of the variables in 
(|14|) . Very accurate results were obtained, depending - obviously - mainly on the 
number of terms included in P(s, t, u): increasing the number of accurately chosen 
polynomial terms resulted in better approximations. The method worked very fine 
for parahelium (and a little less fine for orthohelium) and, after several trials with 
different monomial terms s t m u n (different values of £, in, n), Hylleraas was able to 
obtain the value of W^ c = 24.460 eV with order 6 calculations [26] . 

Soon after the appearance of [26], H.A. Bethe compared this result with that 
obtained with the Hartree method, showing that the overall charge distribution of 
the two electrons given by the Hartree self-consistent field method (with correlation 
effects not included) agreed very well with the predictions of the variational method 

Later in 1929, Hylleraas extended [28] his previous calculations [24], [26] to 
helium-like ions Li + and Bc ++ : 

A most exciting cooperation now started between [me and the fa- 
mous spectroscopist Bcngt Edlen]. No sooner had I found the en- 
ergy of the lithium ion lying in between his limits of error than the 
doubly ionized beryllium ion was on the way, and so it continued 
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with boron 3 plus and carbon 4 plus with me always behind. Sus- 
pecting there were no limits to Edlen's power of producing ions I 
decided to overtake him with a good safety margin by introducing 
even an infinite nuclear charge. This was the origin of the energy 
formula E = -2Z 2 + \Z - e 2 + e 3 /Z - e 4 /Z 2 + ... as counted in 
Rydberg units [23] . 

Hylleraas considered the term \jZr\i in the potential energy as a perturbation 
function, and the factor 1/Z as a variational parameter. The energy eigenvalue E 
and eigenfunction tp were then assumed to be of the forms: 

E = E a Z 2 +ExZ + E 2 +E 3 ^ + (16) 
4> = 4>o + A^, (17) 

and the following result for the ionization potential as function of Z was obtained: 

W Z = Z 2 - ^ + 0.31455-^^. (18) 
4 Zj 

It is interesting to compare the fully numerical evaluation leading to the last two 
terms in Eq. (|18|) with the corresponding analytic (and less accurate) evaluation 
in (JTU]), accounting for a constant 25/128 ~ 0.1953 term. For the helium ionization 
energy there is no particular improvement carried on by this method with respect 
to that in [26], but here Hylleraas pushed the calculations till to order 10, with 
a prediction of W^ c = 24.469 eV, to be compared with the Lyman experimental 
result of 24.467 eV [25]. 

Given the success of the numerical calculations performed by Hylleraas, a more 
mathematical insight into the Ritz variational method, as modified by Kellncr and, 
later, by Hylleraas with the introduction of screening constants, was in order. First 
of all, in 1930, G. Brcit showed that, due to the spherical symmetry of the field, 
the reduction of the general 6-dimensional problem to a reduced 3-dimcnsional one 
could always be performed for S'-states. Then, few months later, C. Eckart proved 
that the quantity E = J ip*Hip&r, where H is the negative energy operator, was 
a lower limit to the term- value of the lowest level of a given spectral series, and 
that the best approximation to that term corresponded to the largest value of the 
integral evaluated over a set of various i\) function. He also showed, however, that 
such approximation was not so good at large distances from the nucleus. 

In the meanwhile, Hylleraas continued to improve his method, and to what 
anticipated in [28] he gave its final form in [31] . Now, the trial wavefunction was 
written as 

tP = C- CT/2 P(cr,T,i/), (19) 

with screened arguments: 

a = 2Zs, t = 2Zt, v = 2Zu. (20) 
Eq. (|18|) was, then, replaced by the more accurate one: 

W,-*-lz + 0*U«-™%» + <»™*, (21) 

with very small numerical coefficients for negative power of Z. In particular, calcu- 
lations of the helium ionization potential were made more precise by three orders 
with respect to the previous estimate, thus obtaining the value Wj — 24.470 eV. 
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The excellent agreement with experiments was very important in putting the 
quantum mechanical approach to atomic spectra, based on many-body Schrodingcr 
equation, on a firm, quantitative ground. As it is apparent from the latest results 
discussed above, the matter rapidly converted from a theoretical problem to a 
numerical one and, in the subsequent years, the increasingly precise studies on two- 
or multi-electron atomic spectra developed into a small industry. The inspection 
of the subsequent editions of the Bethe and Salpeter book [TU], from 1933 to 2008, 
may give a good idea of the evolution of the subject. 

3. MAJORANA CONTRIBUTIONS: THEORY AND NUMBERS 

Majorana was involved in helium in at least two occasions, namely for his second 
[32j and third |33j paper published in 1931, but both papers discuss subsequent 
applications of the helium atom, not directly related to what considered here. Nev- 
ertheless, they also contains key summary remarks worth to be reported. 

Following the seminal work by W. Hcitler and F. London [34] on the formation 
of the H2 molecule, Majorana [32] considered as well the problem of the chemical 
bond, applied to the study of the more intriguing case of the molecular ion Hej" 
[35j . Such a paper gave the first quantitative description of this molecular ion within 
the framework of the Hcitler and London method (and whose results agreed with 
the available experimental data), but, for what concerns us here, it is somewhat 
interesting for the summary made by the author about the wavefunctions for the 
helium atom, to be employed in calculations of observable quantities. 

The eigenfunction of the neutral atom of helium in its ground state 
has been calculated numerically with great accuracy but does not 
have a simple analytical expression [21] [26]. Therefore we need 
to use rather simplified unperturbed cigenfunctions. For example, 
we could assume, as commonly done, for the helium ground state 
the product of two hydrogen-like cigenfunctions with an effective 
Z equal to 1.6 -f- 1.7, depending on the criterium used for the eval- 
uation; if we want to optimize (minimize) the average energy, we 
must then set Z = 2 — 5/16 = 1.6875; if instead we want that 
the diamagnetic constant be in agreement with the experimental 
value and at the same time with the value provided by accurate 
theoretical calculations [2T], we must set Z = 1.60 [52] , 
It was customary for Majorana to relate his theoretical calculations to experimental 
observations [36j . and then to adopt appropriate approximations for the former in 
order to fit the accuracy of the latter, so that such a summary well depict the sit- 
uation around 1930. However, this also explains his reiterated use of hydrogenoid 
wavefunctions and, in particular, his constant search for their possible generaliza- 
tions giving easy but physically meaningful expressions (see below). 

A second paper, published in the same year [33J (but anticipating of few years 
what later discovered by other authors [37]), dealt with the calculation of certain 
doubly excited levels of helium within a variational pcrturbativc approach. Here 
Majorana gave a further generalization of his (unpublished) formula reported in 
Eq. pop , also somewhat reminiscent of some empirical relations found earlier and 
discussed below (see Sect. 13. 2p . The reasoning was as follows. In the zeroth 
approximation (with respect to the mutual interaction between the two electrons), 
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the ground state energy (in Rydberg) is given by W = — 2Z 2 , while in the first 
approximation it will be W = —2Z 2 + 4aZ since the interaction increases as Z. 
Now, 

the second approximation can be evaluated with the method of 
the variation of the unit of length 38J1 that is here equivalent to 
assuming hydrogen- like eigenfunctions with an effective Z* [33J . 

The "more correct expression" he got with this method was the following: 

W = -2Z 2 + 4aZ - 2a 2 , (22) 

the screened charge being Z* = Z — a, that is, just as in Eq. (|T0l) but with 
arbitrary aQ An explicit comparison of what Majorana obtained with the method 
mentioned with the available results was performed too, thus showing the goodness 
of his estimate: 

With this method we find that the value of the ground state is 
— iy/Rh = 729/128 = 5.695, whereas the empirical value and the 
theoretical one obtained by Hyllcraas is — W/Rh = 5.807, the dif- 
ference being less than 2% [33] . 

The work peformed by Majorana behind the results he presented in his papers 
(and not published), however, is far from obvious. Indeed, in his unpublished 
personal study [39] and research [12] notes we find several interesting results and 
methods directly related to the two-electron problem, dating back to 1928-9, some 
of which are apparently preliminary studies for his published papers (especially Ref. 
[33j). In the following we give an account of what obtained by Majorana in those 
years, by sketching some key points. 

Although we do not know precisely when and how Majorana approached the two- 
electron problem (with particular reference to helium), it is a matter of fact that 
in his personal notes we find both purely theoretical contributions and numerical 
calculations, including empirical relations as well. The unpublished material we 
present below does not necessarily follow closely the time evolution of Majorana's 
reasoning. 



3.1. A general variant of the variational method. An interesting modification 
of the variational method was the following [40]. Majorana considered an arbitrary 
function ip in terms of which the wavefunction of the ground state was written as 

ip = aQip + a%H(p, (23) 

H being the Hamiltonian operator. The two variational parameters ao,ai have to 
be determined (for a given tp) by minimizing the energy functional in ([7]), which 
assumes the form 

w _ a o A i + ^QQi^ + ajA 3 
a 2 + 2a a 1 A 1 + a\A 2 ' 



"^Majorana refers simply to the Hylleraas method (see Sect. I2-4Ii : the paper by Fock quoted here 
deals with just a general subject, i.e. the virial theorem in the framework of quantum mechanics, 
where the Hylleraas method is mentioned as an example. 

"^Notc that in Ref. 1331 . Majorana applied this same formula to the ground state as well as to 
other levels (with different values of a). 
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with 



A n — 



(p*H n <pdT 
ip*ipdr 



(25) 



(n = 1,2,3). By introducing the function f(ao,a\) = a^Ai + 2a,QaiA 2 + a\A 3 
WaQ + 2Wa a 1 Ai + WalA 2 , from the minimizing conditions df/da = 0, df/dai 
the set of equations for deducing ao, a\ is obtained: 

a {At -W)+ ai {A 2 - WAi) = 0, 
a (A 2 - WA^ + ai (A 3 - WA 2 ) = 0. 

Now, this is a homogeneous set of equation which admits non trivial solutions for 
a , ai only by requiring the matrix of coefficients to have a vanishing determinant: 

Ai-W A 2 - AiW 
A 2 - A X W A 3 - A 2 W 



= 0. 



(26) 



(27) 



This last condition allows to determine the energy W as the smallest root of the 
corresponding equation, that is 



W 



A 3 - A X A 2 + ^(A 3 - A 1 A 2 y - 4(^3 - A\){A 2 - A\) 



(28) 



2(^2-^) 

In the Quaderni Majorana did not apply such a method to come up with numerical 
predictions about the ground state energy of helium, but it is nevertheless inter- 
esting to compare such predictions with already known results. For example, by 
choosing the form of tp as that in <j8j) , we have obtained Wf~ c = 21.62 eV (for 
k = 1.6711), which is certainly not a good result when compared already with what 
obtained by Majorana himself with the standard variational method for the same 
trial wavefunction (see Sect. |2~2]) . i.e. Wf c = 22.95 eV. 

This is probably the reason why Majorana further generalized his method, by 
writing the wavefunction: 

V' = a ip + a x H<f> + a 2 H 2 ip + . . . + a n H n Lp (29) 

for a given (p and a general index n, for which the energy to be minimized is (with 
the same notations as above): 



^ Oja^ A i+ k+i 



W 



i,k=0 



(30) 



i 

l W will be the smallest root of the following equation:" 



E 

i,k=0 



A 1 


- W 


A 2 


- A X W 


A n - A n -iW 




A 2 


- A X W 


A 3 


- A 2 W 


A n+ l ~ A n W 




A 3 


-A 2 W 


A 4 


- A 3 W 


A n +2 - A n+1 W 


= 


A„ 




A n 


H-A n W .. 


A 2n -i - A 2n - 2 W 





(31) 



Obviously, by increasing n, the calculations give better and better results but, as 
well, become more and more difficult. However, the point on which Majorana 
focused was not that but, rather, the emergence of divergencies: 
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Often, this procedure does not converge, because, starting from 
a given value of n, quantity H n tp exhibits too many singularities, 
which forces us to consider only combinations of the form 

ip = a ifi + aiHip + . . . + a n -\H n ~ x <f. (32) 

The inclusion of additional terms is not useful, since the corre- 
sponding a coefficients would necessarily vanish. 

He referred to the calculations of integrals for r± — 0, r2 — or r\i = 0, whose 
degree of divergency increases with incrasing n in H n , given the form of the Hamil- 
tonian operator. It is quite interesting the connection made by Majorana between 
the form (for a given ip) of ip in ([32]) , dictated by choosing a suitable n, and the 
divergency properties of the system, as described by its HamiltonianQ 

3.2. Empirical relations. Seemingly, the first numerical studies on the two-electron 
problem performed by Majorana were aimed to find empirical relations for the 
ground state energy [41] . He was evidently aware of the difficulties in obtaining 
theoretical predictions^ so that the first numerical efforts were likely devoted to 
find suitable empirical expressions for W . According to Majorana himself, several 
of these relations lead to unsatisfactory results, but it is nevertheless intriguing the 
reasoning behind themQ 

Let us consider a two-electron atom with charge Z in its ground 
state. We denote by a = < 1/n > = < l/r 2 > the mean value of 
the inverse of the distance of each electron from the nucleus, and 
with b =< l//'i2 > the mean value of the inverse of the distance 
between the two electrons. Expressing the distances in electronic 
units and the energy in Ry, we have 

W = -2aZ + b, (33) 

since the energy is equal to half the mean value of the potential 
energy. If we now consider an atom with atomic number Z + dZ, 
perturbation theory gives 

dW = -4adZ, (34) 

and thus we have two equations in the three unknown Z functions 
E,a,b HE]. 

Interestingly enough, Majorana's strategy for dealing with generic two-electron 
atoms with atomic number Z was to consider this quantity as a continuous variable, 
a small change of which induces a change in the ground state energy, the latter one 
being evaluable with the aid of the perturbation theory. Now, the problem was to 
have two equations with three unknowns, so that the necessary additional relation 
was introduced as an empirical formula: 



^Recall that such contribution dates around 1928-9. 

^Apart from the references in I32| and 1331 . we do not know precisely what papers were read 
and studied by Majorana, but from the arguments and, especially, from the symbols he used in 
his notes 1391 1121 , we can safely deduce that he was well aware of a large part of the literature 
quoted above. 

^In what follows, we use consecutive labels for the equations in the quotations from the original 
Majorana's manuscript; we also conform, when necessary, the original notation to the present one. 
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We now add another empirical relation between a and b, which is 
presumably a good approximation: 



b = (2Z - 2a) (2a - Z). 



(35) 



This relation can be deduced from the following considerations. For 
sufficiently high values of Z ', perturbation theory gives 




(36) 



but, on the other hand, 



b = lz + 



(37) 



so that, from Eq. (l33l) . 




(38) 



which, in first approximation, satisfies Eq. Q35p. For very small 
values of Z we can consider that the first electron is next to the 
nucleus, while the other one is practically at an infinite distance; 
then we have a ~ Z/2, b ~ 0, and Eq. ([33]) is again satisfied. We 
finally assume that it is also a good approximation for intermediate 
values of Z. 

The next step was, then, simply to substitute Eq. ([33]) into Eqs. (f3"3"| and (f3~4"|). 
thus obtaining dZ = da, and "since we know the value of a for infinite Z, we deduce 
a = Z — 5/16." The ground state energy is, then, accordingly, 



to be compared with the perturbative result in Eq. © and subsequent improve- 
ments. Majorana noticed that such formula could be used only for Z > 5/8, while 
for Z = 5/8 the quantity b in (|35|) vanishes. Moreover, "the procedure used here 
is not very satisfactory, since for very small Z the quantity b would vanish faster 
than a first-order term and it would become negative." A comparison with data for 
the helium atom revealed that the predicted value for W was in excess of 1.13 eV 
with respect to the experimental value (25.59 eV instead of 24.46 eV), but here the 
intriguing thing is that Majorana applied the above formula also to hydrogen: "For 
the hydrogen atom [Z = 1) we find instead W = —1.141, from which the ionization 
potential would be 1.91 (electron affinity). "0 

Given the poor agreement with the experimental values, Majorana also explored 
other relations for b as functions of Z, such as, for instance, 



W = -2Z 2 + -Z 



4 



25 
64' 



(39) 




(40) 



or 




(41) 



or even 



b = 



-Ze~ k / Z 



(42) 



"Note that the predicted electron affinity of the hydrogen atom, i.e. the difference between 
the ground state energies of the neutral atom and the once-ionized atom, is more than twice the 
actual value. 
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where k is a parameter to be determined under suitable conditions. For example, 
from the last relation, he obtained the following expression for a: 

where it is evident the contribution from any Z ranging from to oo. 

Since for small Z we must have a ~ Z/2, wc can choose k such that 



f-5Z + 5fc v 1 
Jo 16^ 3 2' 



i.e., k = 5/4. However, in this way we obtain a bad approximation. 
Indeed, we would have L.l 



(44) 



E = -& - Z*e- 1M ' Z , (45) 

and for helium (Z = 2) we would get E = —5.14, which is a value 
far from the experimental one. 

Majorana was well aware that the introduction of such empirical relations was 
completely arbitrary ( "there is no a priori reason to prefer one or the other" ) , so 
that he soon abandoned these calculations. Here, however, it remains the basic, 
novel assumption that Z could be treated - in an "effective" way - as a continuous 
variable, and that contributions to the ground state energy of a given atom (with 
given atomic number) come also from any other Z . 



3.3. Helium wavefunctions and broad range estimates. The conclusions 
achieved by Slater in 1928 [OEj about the form of the unperturbed wavefunction 
and the role played by the operator H in the perturbative method (see Sect. 12. ip . 
were analyzed repeatedly by Majorana [42] , who considered as well several different 
forms of the helium wavefunction as generalization of (|4j). along the same line con- 
sidered by Slater (who introduced, in particular, the simple form in ((6])). However, 
from the given wavefunctions, Majorana was able to deduce also (broad) ranges 
within which the value of the helium ground state energy would have lain, thus 
quantifying in a definite manner what generally deduced by Slater. 

Majorana realized very early what later employed extensively by Hylleraas (but 
already latent in Slater's calculations), namely that the helium wavefunction de- 
pends only on three metric quantities r\,T2,ri% (sec also Sect. I3.4j) . Probably 
inspired by Slater's Eq. ([TB")) , Majorana firstly introduced the following general 
form: 

if> = e~*, 

2ri + 2r 2 - -7*12 + a{r\ + r\) + br\ri + cr\ 2 + d(n + r 2 )ri 2 (46) 
P 1 + e(ri + r 2 ) + fr 12 

where a, b, . . . , / are numerical coefficients. The energy eigenvalue equation ([1} was 
then written (in suitable units) as 



Lip = Xip 



(47) 
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with 



4 4 2 



2 



L = — + + V 

n r 2 ri 2 

4 4 2 d 2 d 2 d 2 2 d 2 d 

2 



n r 2 r 12 dr\ dr\ dr\ 2 n dri r 2 dr 2 r 12 dr 12 

d 2 d 2 

+2cosai-- — h2cosa 2 -^ — r — , (48) 

aridr 12 dr 2 dr 12 

where «i[«2] is the angle between ri[r 2 ] and r i2 . However, he soon realized that 
calculations would become quite difficult with such a general form, so that he shifted 
his attention to the simpler form 

V^o =e~ 2ri - 2r2+ ^ ri2 . (49) 

Though, it should be noted that such a choice was not accidental (or dictated only 
by computational simplicity): indeed, he deduced the relations to be hold among 
the a — f coefficients in order to have from Eq. (|4"f)|) the same result (see below) 
as coming from Eq. (|49|) 1^ By substituting (|49l) into the eigenvalue equation (l47l) . 
Majorana obtained simply: 

17 

A= — — 2 cosai — 2 cosa2. (50) 

From this expression, by allowing the cosines to take any value in their domain, he 
deduced the range inside which the (negative) energy eigenvalue has to lie, namely 
4.5 < -VK/Rh < 8.5. 

When compared with the value of 5.807, this is indeed a very broad range, but, 
as said above, it should give a quantitative idea of the results obtainable within 
the perturbativc method, as induced by the particular form of the Hamiltonian 
operator on different unperturbed wavef unctions. Therefore, Majorana passed to 
consider an "approximated" wavefunction, that is 

^={l + \ri^je- 2 ^- 2r \ (51) 

Note that such expression approximates that in (|4^|) (and, then, that in (HBJ)) f° r 
small r 12 , i.e. when the two electrons are close each other: it can be then regarded 
as a correction to the simple one-electron case. By repeating the same calculations 
as above, Majorana obtained 

1 2 

A = 8 j ; (cos ot\ + cos a 2 ) , (52) 

1 + 57*12 1 + |ri2 

and thus 3 < — M^/Rh < 8. The comparison of this range with that obtained 
previously already quantifies to a certain extent the result obtained by Slater [T3] 
and addressed above, but Majorana did not limit himself to just this two forms of 
the possible wavef unctions. Without entering into computational details [12], we 
only mention a different approximation of (|46[) , namely 

ip = 1 - 2ri - 2r 2 + -r u + a(r\ + r\) + br\ 2 + cr x r 2 + d{n + r 2 )r 12 + ... (53) 



"For technical details, see 1421 . 
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for small T\,T2, r\i and a generalization of ([5T 



l + \ Tl2 



-(2-2a)r 2 -(2-2a)n-2r 2 



1 + 2or 2 1 + 2ott 



(54) 



(obviously symmetric under the exchange 1 f-> 2), with the parameter a parame- 
terizing an effective nuclear charge different for the two electrons, probably due to 
the different position of the two electrons considered by Majorana with respect to 
the nucleus. 



3.4. A simpler alternative to Hylleraas' method. As we have in Sect. 12 .41 a 
sensible improvement in the predictions obtained with the variational method came 
with the recognition by Hylleraas that the wavefuntion of the two-electron problem 
depends only on three coordinates (instead of six), namely the magnitude of the 
two vectors ri,r2 and the angle 9 between them or, rather, ri,r2,ri2. In addition 
to the perturbative calculations just discussed, and likely before the appearance of 
the Hylleraas papers [24j [26] (or, in any case, independently of them), Majorana 
implemented such a result in one of his (unpublished) variational studies about the 
ground state energy of two-electron atoms |43j , aimed at obtaining an improvement 
of his general result in Eq. (TTUl) . Instead of the wavcfunction in © , he considered 
the function 

?P = e- kri e- kr2 e iri2 (55) 

with arbitrary k and £ to be determined: "we will certainly obtain a better approx- 
imation." The energy in (|7|) was, then, evaluated by noting that: 

Hip = -2k 2 ^ - 2£ 2 tp + 2k£ipai + 2k£^a 2 - 2^—^ - 2^— -tp + 2 - 4 % , (56) 

n r 2 ri2 

where a\ [a 2 ] is the cosine of the angle between ri [r 2 ] and ri 2 . However, as 
subsequently done also by Hylleraas, the angular quantities were replaced by only 
metric ones by means of the formulae: 

r\ + r\ 2 - r| r| + r\ 2 - r\ 

ai = , a 2 = . (57 

2riri2 2r 2 ri2 

For some reason, apparently Majorana did not completed the full calculations (see 
|43| ) . which would have resulted in the following expression for the energy to be 
minimized: 

_ 2(k - I) [8fc 3 + fc 2 (5 - 11 - 16Z) + 4kl(Z + £ - 1) - £ 2 (£ - I)] 
W ~ 8 fc 2 - 5k£ + £ 2 ' (58) 

As expected, the minimization procedure gives a better resulQ f° r the helium 
ground state energy with respect to the £ = case. However, despite the apparent 
complicated formula in (|58p . the striking fact is the accuracy of the Majorana 
result, obtained just at first order, Wf c = 24.09 cV (corresponding to k — 1.8581, 
£ = 0.2547), when compared with the analogous one produced by Kcllner at fourth 
order, Wf~ c = 23.75 cV, and even with that obtained by Hylleraas with calculations 
at order 11, Wf c = 24.35 eV. 



^In particular, Majorana deduced the numerical values for the coefficients by requiring that 
ij> and its first derivative have a node at the same position when the two-electron system collapses 
into a one-electron one, i.e. n = (or = 0) and r\i = 0. 

n -iy/Rh = 5.779 instead of -W/Rh = 729/128 = 5.695. 
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4. Conclusions 

The early history of quantum mechanics was indissolubly related to the successful 
description of the helium atom, with particular reference to the prediction of its 
ground state energy. The basic difference with the description of the simplest atomic 
system, the hydrogen (or hydrogenoid) atom, is that the two-electron problem does 
not admit an analytic solution, so that quantitative predictions about it are de 
facto mediated by the approximation method employed to get numerical estimates, 
to be then compared with the experimental observations. 

As seen in Sect. [2j though depending on the chosen unperturbed system, per- 
turbative calculations performed since 1926 did not give very accurate predictions, 
even by increasing the perturbativc order, the reason for this (envisaged by Slater) 
being due to the form of the Hamiltonian operator for the system at hand. Such an 
occurrence led physicists to adopt different approaches, the first of which was the 
variational method early introduced by Ritz. Calculations performed by Kcllner 
in 1927, indeed, did improve perturbative results, but the interesting novel idea 
introduced by him was that of an effective nuclear charge experienced by one of 
the two electrons in the field of the nucleus plus the other electron, which shields 
to some extent the nuclear charge. 

In the meanwhile, at the end of 1927 Hartree elaborated his self-consistent field 
method aimed at describing many-eletron atoms, whose basic assumption was just 
that such atoms could be described by wavefunctions given as the product of one- 
electron wavefunctions of individual electrons in the field of a corrected ("self- 
consistent") central potential. Although criticized by many, mainly due to its 
"unusual" line of reasoning, the Hartree method gave a prediction for the helium 
ionization potential very close to the experimental determinations. 

Important novelties on the subject came early in 1928 by means of Hylleraas, who 
realized that the absolute orientation of the atom in the space was irrelevant, this 
leading to an enormous yet unexpected simplification into the variational problem, 
described only in terms of three variables defining the shape of the electron-nucleus 
triangle, instead of the six spatial coordinates of the two electrons. With this sim- 
plification, Hylleraas was able to perform order 11 calculations (instead of only 
order 4 ones by Kcllner), which provided an excellent approximation for the ion- 
ization energy of helium. The method and reasoning elaborated by Hylleraas, with 
the subsequent introduction of what will become known as Mandelstam variables 
s,t,u, were worth to be easily generalized, and the later introduction of screening 
constants further improved the predictions. The quantum mechanical two-electron 
problem converted, around 1930, just into a numerical one, whose accuracy de- 
pended only on the degree of approximation used. 

This known history of the solutions given to the two-electron problem has been 
here supplemented by several, unpublished results contributed by the Italian physi- 
cist Ettore Majorana around the same years, as described in detail mainly in Sect. 
[3] Indeed, almost all of these results did not convey (even indirectly) in his two 
papers on the topics, published in 1931, so that they remained unknown until now. 
Due to their simple deduction, Majorana's reproduction of standard perturbative 
and variational results, but generalized to arbitrary Z, has been reported as well 
(in Sect. [5]) for pedagogical reasons: they are summarized by the formulae in Eqs. 
(0 and (fTO)) . respectively. 
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From a purely theoretical point of view, the most intriguing contribution is a 
variant of the variational method, developed by Majorana in order to take directly 
into account, already in the trial wavcfunction, the action of the full Hamiltonian 
operator of the given system. From the general wavcfunction in Eq. (|29p he was, 
then, able to obtain a detcrminantal equation (|31l) from which the ground state 
energy could be deduced. Though this method was originally devised for helium 
(or, more generally, for the two-electron problem), it is well suited for application 
to any variational problem and, as such, it is worth to be applied in different fields 
of present day research. 

Instead, coming back to the problem of the ground state energy of two-electron 
atoms, it is particularly interesting to mention an application of the perturbative 
method where the atomic number Z was treated effectively as a continuous vari- 
able, contributions to the ground state energy of an atom with given Z = Zq 
coming also from any other Z ^ Zq. Its application to several empirical relations 
aimed at finding accurate predictions - without the use of complicated unperturbed 
wavefunctions - has been described in detail in Sect. 13.21 

The use of particular wavefunctions, along a line of reasoning similar to that 
adopted by Slater in 1928, was as well considered by Majorana. In addition to the 
conclusions achieved by Slater, regarding the accuracy of the perturbative method 
(obtained, however, in a more transparent and quantifiable manner, with the de- 
duction of broad range estimates), he introduced the concept of an effective nuclear 
charge different for the two electrons, thus generalizing previous results (see Sect. 

Finally, we have discussed a simple alternative to Hylleraas' method which, as 
explained in Sect. 13.41 led Majorana to obtain first order numerical results compa- 
rable in accuracy with those by Hylleraas at order 11 (and certainly better results 
with respect to those by Kellner at fourth order), just by choosing a suitable (and 
quite "obvious" ) trial wavcfunction. 

Summing up, if the role played by Majorana's contributions in the known liter- 
ature on the two-electron problem was, in a sense, certainly subdominant, never- 
theless his unpublished results add new non negligible pieces to the history of early 
quantum mechanics. But, especially, their pedagogical relevance in undergraduate 
and postgraduate courses, as well as the importance of novel methods in current 
researches in quantum physics, as described in the present paper, is yet to be fully 
discovered and exploited. 
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